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We show that under the influence of pure vacuum noise two entangled qubits become completely 
disentangled in a finite time, and in a specific example we find the time to be given by In ^^-^^ 
times the usual spontaneous lifetime. 
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Superposition and entanglement are two basic features 
that distinguish the quantum world from the classical 
world. While quantum coherence is recognized as a major 
resource, decoherence due to the interaction with an envi- 
ronment is a crucial issue that is of fundamental interest 
0, 0, 0, • When coherence exists among several dis- 
tinct quantum subsystems the issue becomes more com- 
plicated because, along with the local coherence of each 
constituent particle, their entanglement brings a special 
kind of distributed or nonlocal coherence. It is this dis- 
tributed coherence that really matters in many impor- 
tant applications of quantum information |^ Con- 
sequently, the fragility of nonlocal quantum coherence 
is recognized as a main obstacle to realizing quantum 
computing and quantum information processing (QIP) 
Apart from the important link to QIP realizations, 
a deeper understanding of entanglement decoherence is 
also expected to lead to new insights into quantum fun- 
damentals, particularly quantum measurement and the 
quantum-classical transition ITol . [Tl| . Although quan- 
tum decoherence has been extensively studied in recent 
years, it remains unclear how a local decoherence rate is 
related to a nonlocal disentanglement rate when a multi- 
particle quantum state is in contact with one or more 
noisy environments. 

Therefore, a deep understanding of the decoherence 
in any viable realization of qubits is desirable and it is 
surprising that few if any fundamental treatments exist 
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FIG. 1: Schematic illustration of a set-up in which two atoms 
are located inside two spatially separated cavities A and B. 
The two atoms are initially entangled but have no direct in- 
teraction afterwards. 



of decoherence that include the dynamics of disentangle- 
ment on better than an empirical or phenomenological 
basis. 

Here we consider two initially entangled qubits and 
examine the dynamics of their disentanglement due to 
spontaneous emission without phenomenological approx- 
imation. There is perhaps no simpler realistic bipartite 
model in which all of the effects of quantum noise can 
be considered fully analytically. We show that decoher- 
ence caused by vacuum fluctuations can affect localized 
and distributed coherences in very different ways. As 
one surprising consequence, we show that spontaneous 
disentanglement may take only a finite time to be com- 
pleted, while local decoherence (the normal single-atom 
transverse and longitudinal decay) takes an infinite time. 

To make our model and results concrete, we restrict 
our attention to two two-level atoms A and B coupled 
individually to two cavities which are initially in their 
vacuum states (see Fig. In the general framework 

of system-plus-environment, the two two-level atoms are 
identified as the system of interest, whereas the two cav- 
ities serve as the environments. The interaction between 
each atom and its environment results in the loss of both 
local coherence and quantum entanglement of the two 
atoms. In its simplest form such a model may be for- 
mulated with the following total Hamiltonian, which is 
given by (we set H = 1): Htot = -^at + ^^int + ^^cav, where 
the Hamiltonian of the two atoms -ffat, the two cavities 
i/cav and the interaction iJjnt are given by 
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where gk, fk are coupling constants and denotes the 
usual diagonal Pauli matrix, and the standard 2-qubit 
product basis is given by: 
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where | ± ±)ab = |±)a ® |±)b denote the eigenstates of 
the product Pauli spin operator cr^ ® erf with eigenvalues 
±1. The total Hamiltonian, given by equations Q-©, 
provides us an important solvable model of the atom- field 
interaction in quantum optics. 

Suppose that initially the atoms are entangled with 
each other but not with the cavities, i.e., we assume that 
at i = the two atoms and the cavities are described by 
the product state. 



Wab(^\0)a\0)i 



where |'0)as is the entangled initial state of the two 
atoms and |0)a|0)s is the vacuum state of two cavities. 
For simplicity, we will not take into account the spatial 
degrees of freedom of the two atoms. The convolution- 
less master equation of the system of two atoms can be 
obtained as follows 12]: 
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where H'^^ (t) is the system Hamiltonian modified to take 
into account the Lamb shifts 
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Kdt) = ^{^A + Fi{t))af + -{UB + Gi{t))af (7) 

and the coefficients Frj and Grj are the real and imag- 
inary parts of F{t) and G{t), which are given by 
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dsa{t — s)h{s) 



dsf3{t ~ s)c(s) 



Note that the cavity correlation functions are 
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iujAbis) + / dXa{s - X)b{X) = (12) 
Jo 

iuBcis) + I dXl3{s ~ A)c(A) = (13) 



and the functions b{s) and c(s) are the fundamental so- 
lutions of the equations of motion: 
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The master equation © is extremely useful for the 
study of decoherence, and for the purpose of disentangle- 
ment analysis it will be very convenient to find an explicit 
expression for its solution. In the interaction picture, 
where 

pit) = e*/o ^^t(s)d«p(t)e-^/o ^atM-i-^ (14) 



the general solutions of equation © can be described in 
terms of a Kraus representation ^ 0, 0| . As can 
be seen below, the Kraus representation allows a very 
elegant analysis of the disentanglement time for arbitrary 
states. Precisely, for any initial state p(0), the density 
operator at t can be expressed as|l3 
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p(0 = ^i^p(t)p(0)i^t(i)^ 



(15) 
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where the Kraus operators K^it) satisfy X^u^u-^'^m = ^ 



fj.\~y J fJ, 

for all t. The Kraus operators for this model are gi" 
by 
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and the time-dependent Kraus matrix elements are 



J Ait) = exp 
7b (t) = exp 
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UAit) = ^jl-llit), ujBit) = ^l-ilit) (22) 

With the preceding discussion, we are now in a posi- 
tion to determine both the local decoherence rate and 
the disentanglement rate. Multiple interpretations of the 
term decoherence in the literature can lead to confusion, 
so we will use global or non-local decoherence (or disen- 
tanglement) here when we refer to loss of bipartite entan- 
glement. The terms local decoherence or local relaxation 
will refer to longitudinal and transverse decay of single- 
atom density matrix elements. In the present example 
local and non-local decoherence both arise from the ef- 
fects of spontaneous emission and in that sense are not 
independent. 

We begin with the coherence decay of a single qubit un- 
der the master equation © . The local decoherence rates 
of the qubits can be estimated from the reduced density 
matrices p"^ = Ttb{p} and p^ = Tia{p}- The lo- 
cal decoherence rates are determined by the well-known 
Bloch equations with general time-dependent functions 
F(s),G(s). For example. 



[Fit) + iua] (a^) 



(23) 



and (ct^) 



(T_)*, where a_ = (j^®F Similar equations 



hold for cr^. Hence we have 



\{aiit))\ ^ |(a^(0))|e-/o'i^^^H(^), 
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Given these equations, local decoherence behaviors are 
determined by the character of the functions Ffi{t) and 
Gn{t) which we always assume to be positive functions 
asymptotically. In the familiar Born-Markov approxima- 
tion one has purely exponential decay with rates given 
by F{t) rA/2 and G{t) Tb/2, where the T's are 
the Einstein A coefficients for the two-level atoms in the 
cavities. 

The comparison of interest is with the disentangle- 
ment rate. Since entanglement decoherence processes 
are most generally associated with mixed states, we will 
use Wootters's concurrence to quantify the degree of en- 
tanglement [l^ . The concurrence is conveniently de- 
fined for both pure and mixed states. Let p be a den- 
sity matrix of the pair of atoms expressed in the stan- 
dard basis 10}. The concurrence may be calculated ex- 
plicitly from the density matrix p for qubits A and B: 
C{p) — max (O, VA7 — — ~ ) j where the 
quantities A.; are the eigenvalues of the matrix C: 



(26) 



arranged in decreasing order. Here p* denotes the com- 
plex conjugation of p in the standard basis Q and ay 
is the usual (pure imaginary) Pauli matrix expressed in 
the same basis. It can be shown that the concurrence 
varies from C = for a disentangled state to C = 1 for 
a maximally entangled state. 

We now show two categories of result for entanglement 
decay. In the first more general category we show that, 
for all entangled (possibly mixed) states, entanglement 
decays not only more rapidly than the fastest decoher- 
ence rate of an individual qubit, but at least as fast as 
the sum of the separate rates. In the second category we 
present a sharper result in a specific mixed state exam- 
ple, in which the entanglement goes exactly to zero in a 
finite time and remains zero. Both categories of result 
are a consequence of normal spontaneous emission. 

For the first category of result, let us note that the 
concurrence C(p) is a convex function of p ,,16j . From 
(|15|l . one immediately has 



C{pit))<J2c{K,piO)Kl)., 



(27) 



where Kf^ are defined in equations H16() to H19|) . Let us 
consider a typical term C{K^p{0)Kj^) in l|77|l and denote 
it by Pfi = Kfj^p{0)Kj^. From the definition of concurrence, 
it can be proved that 

C(pi) = e-/o(-^«(^)+^«(^»''^C(p(0)), (28) 
C(p^) = 0, Ai = 2,3,4. (29) 

Then the inequality H27() immediately leads to, 

Cip{t)) < e--fo(-P«(-)+G«(«))'i-c(p(0)), (30) 

which establishes the first result mentioned. It is not 
difficult to show that the upper bound is the minimal 



upper bound. To treat the disentanglement process in 
general and more completely requires a discussion of the 
asymptotic behavior of the functions Fii{s) and Gr{s), 
which is beyond the scope of the present paper. 

In what follows we develop the second result mentioned 
above. We show that within the general result there are 
very unusual and striking specific consequences. One ex- 
ample shows that, within the general exponential char- 
acter evident in H3(J|I . disentanglement can be completed 
in a finite time while the local decoherences need an infi- 
nite time. Let us assume that the initial density matrix 
is only partially coherent, but include an arbitrary de- 
gree of non-local coherence of a familiar type (one of the 
atoms is excited, but it is not certain which one). This 
is easily expressed in the following form 



/a{t) 

b{t) z{t) 

z*{t) c{t) 

V d(i)/ 
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where the factor 1/3 is for notational convenience. The 
concurrence for this density matrix is C = | max{0, \z\ — 
Vad}. For simplicity, we consider an important class of 
mixed states with a single parameter a satisfying initially 
a > 0, d = 1 — a, and 6 = c = 2; = l, so then initially 
C(p(0)) = |[1 - y/a{l - a)]. For t > the matrix ele- 
ments are given by 
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= lA + lA^la, 
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= 1b+ ^Alsa, 
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= 1 — a + iLj\ + LJ% + uj^LU^a, 
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To simplify the calculations, we use the Markov limit 
results and assume the cavities are similar so 7a = 7b = 
7 = exp[— rt/2], and loa — = lo = \/\ — exp[— Ft]. 
The concurrence for the density matrix p(t) is given by 



C(p(<)) = -max{0, 7^/(0} 



(37) 



where f{t) — 1 — y/a{l ~ a + 2lj^ + uj'^a). The sufficient 
condition for the concurrence 137|l to be zero is 



l~a{l-a + 2LU^ + Lu^a) < 0. 



(38) 



The simplest case is a = 1 and from this we can easily 
show the surprising result that the density matrix l|31|) 
has a finite disentanglement time. That is, C{p{t)) = 
for all t > td, where td is very finite: 



td 



2-1- V2 



(39) 



In response to this surprising result, a natural ques- 
tion will be: does spontaneous emission cause all initially 
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FIG. 2: This figure sliows the entanglement decay via spon- 
taneous emission starting from the initially entangled states 
1311 1 with a between zero and 1, d = 1 — a and b — c = z = 1. 
Finite-time complete disentanglement takes place for a > 1/3, 
and when a < 1/3 disentanglement of the initial state is only 
completed asymptotically. 

entangled two-qubit states to disentangle at some finite 
critical time? The answer is no. To see this, we consider 
the entire range of different a values, and plot the con- 
currence decay in Fig. |21 The figure shows that for all a 
values between 1/3 and 1 concurrence decay is completed 
in a finite time, but for smaller a's the time for complete 
decay is infinite. The different behaviors exhibited over 
the allowed range of a values in Fig. |2]show that our two- 
level atom model behaves qualitatively differently from 



the continuous variable two-atom model discussed in 
by Dodd and Halliwell. In that case all initially entan- 
gled states become separable after a finite time (see also 

E3). 

To summarize, we have shown for the physically funda- 
mental and unavoidable process of spontaneous emission 
that nonlocal disentanglement times are shorter than lo- 
cal decoherence times for arbitrary entangled states (pure 
or mixed). We based our results on perhaps the sim- 
plest realistic decoherence scenario in which two entan- 
gled qubits individually interact with vacuum noise. The 
model allows an exact analysis and also shows, remark- 
ably, that complete disentanglement can be reached after 
only a finite time, whereas more familiar local decoher- 
ence processes take an infinite time to be complete. We 
believe our results are of generic nature. Undoubtedly a 
deep understanding of the relation between decoherence 
and disentanglement will be of importance for both the 
foundation of quantum mechanics and practical quantum 
information applications. 
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